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Abstract 

We consider the problem of a single level quantum dot coupled to the edge of a one-dimensional Luttinger liquid wire by both 
a hopping term and electron-electron interactions. Using bosonization and Coulomb gas mapping of the Anderson- Yuval 
type we show that thermodynamic properties of the level, in particular, its occupation, depend on the various interactions 
in the system only through a single quantity — the corresponding Fermi edge singularity exponent. However, dynamical 
properties, such as the level density of states, depend in a different way on each type of interaction. Hence, we can construct 
different models, with and without interactions in the wire, with equal Fermi edge singularity exponents, which have identical 
population curves, although they originate from very different level densities of states. The latter may either be regular 
or show a power-law suppression or enhancement at the Fermi energy. These predictions are verified to a high degree of 
accuracy using the density matrix renormalization group algorithm to calculate the dot occupation, and classical Monte 
Carlo simulations on the corresponding Coulomb gas model to extract the level density of states. 
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1. Introduction 

The behavior of low-dimensional electronic systems 
has been in the focus of many experimental and the- 
oretical studies in recent years. Such systems are im- 
portant both because of the fundamental interest in 
the strongly correlated physics they exhibit, as well as 
their role as the building blocks for creating nano-scale 
devices. The one-dimensional case is a particularly in- 
teresting one. When no symmetry is spontaneously 
broken and the one-dimensional system is metallic, its 
low energy dynamics is described by the Luttinger liq- 
uid (LL) theory [1]. The latter offers one of the clear- 
est realizations of non Fermi liquid physics. These sys- 



tems have been experimentally realized in a variety of 
ways, including narrow quantum wires in semiconduct- 
ing heterostructures, metallic nanowires, and carbon 
nanotubes. A closely related concept is that of chiral 
LLs, which describe the physics of edge states in the 
fractional quantum Hall effect [2]. A natural question 
which arises is the effect of impurities on these sys- 
tems, either naturally occurring or artificially intro- 
duced (e.g., quantum dots and anti-dots). While this 
topic has been investigated for some time, most of the 
previous works are restricted to the study of transport 
properties [1-5] , while other phenomena have been only 
occasionally considered [6-13]. 

In this paper we investigate one of the simplest sys- 
tems of this kind, namely a single level attached to the 



Preprint submitted to Physica E 



20 January 2011 



end of a LL (or, equivalently, a level in the vicinity 
of a chiral LL). We include short-range interactions 
between the charge of the level and the charges in its 
neighborhood in the wire. Although transport prop- 
erties are, of course, not relevant here, many other 
phenomena can be investigated. In this paper we 
compare thermodynamic properties, in particular, the 
level population (which can be measured by placing a 
quantum point contact in the vicinity of the level), to 
dynamic properties, such as the level density of states 
(LDoS; this can be probed by tunnel spectroscopy). In 
section 2 we demonstrate, using the Anderson- Yuval 
approach, that while thermodynamic properties are 
universal, and are affected by the various interactions 
only through a single parameter (identified as the 
Fermi edge singularity exponent), dynamic properties 
are sensitive to the specific physics of the different 
interaction types. Based on this analysis, in section 
3 we construct different systems, some LLs and some 
not, which are tested numerically to have the same 
dependence of the population on the level energy, al- 
though their LDoSs (of which the population is an 
integral) are very different. These results are discussed 
and summarized in section 4. 



2. Model and Coulomb gas analysis 

The system described in the Introduction can be 
modeled by the Hamiltonian H — H„ + Hi + //i w - 
The first term is the wire Hamiltonian, described by 
the standard Tomonaga-Luttinger model. In bosonized 
form, it is given by [1]: 



oc 

Hw = ^I {^[V0(z)] 2 + s[v<Kx)] 2 }dz, 



(1) 



where the bosonic fields 6(x) and 4>{x) obey the com- 
mutation relation [6(x), <f>(y)] = i7r6(a; — y), and the 
boundary condition 6(0) = 0. Here g and v are the 
usual interaction parameter and excitation velocity, re- 
spectively, and Q(x) is Heaviside's step function. The 
second term of the full Hamiltonian describes the level: 
Hi — £(,d)d, where (ft (d) is the level creation (annihi- 
lation) operator and £o is its energy. Finally, the level 
and the wire are coupled by: 



fliw = - [7T w dty(0) + H.c] + 



(2) 



The first part of H\ w describes the level-wire hopping, 
parametrized by a tunneling matrix element 7T W , while 
the second part, in which the colons denote normal 
ordering, is a local level-wire interaction of strength 
Wiw The electronic annihilation operator at the wire's 



edge is given by 4>(0) = xc 



<K0) , 



/27m, using Majorana 



Fermi operators x an d a short distance cutoff (e.g., a 
lattice spacing) a. 

Following Anderson and Yuval's method [14], any 
quantity of interest is expanded to all orders in 7T W - 
This results in a series of correlation functions, which 
are calculated for 7T W = [15,16]. Because of the level- 
wire interaction, there is a potential at the edge of the 
wire which flips between ±Wi w when each hopping oc- 
curs, i.e., we have a sequence of Fermi edge singularity 
events [17]. 

Let us first discuss the thermodynamic properties of 
the model, e.g.: the level population (to be denoted by 
Tiievei), entropy and the specific heat. These can be ex- 
pressed through the partition function Z of the model 
and its derivative with respect to the parameters of the 
system, such as the level energy £o and the tempera- 
ture T. In the Anderson- Yuval approach, the partition 
function acquires the form of a grand canonical parti- 
tion function of a classical system of particels (hopping 
events) residing on the imaginary time axis of the orig- 
inal quantum model (i.e., on a circle whose circumfer- 
ence is the inverse of the temperature T of the level- 
wire system) . We assign to a particle corresponding to 
hopping of an electron from the level to the wire a pos- 
itive charge, and to a particle describing the reverse 
process a negative charge. Thus, the charges must be 
alternating in sign, and their total number has to be 
even. Denoting the position of the i'th particle by n 
and the sign of the charge of the first one by s, we ob- 
tain: 
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so that the particles have a fugacity yjTo^o/ir, where 
£o is a short time (ultraviolet) cutoff and To is the level 
width (an expression for which is given below). The 
Hamiltonian of the classical system reads: 

H CG ({n},s) = 



(4) 



The first part of this expression describes an interac- 
tion between the particles, which is similar in form 
to the electrostatic interaction between charged rods, 
giving Eqs. (3-4) the name "Coulomb gas expansion". 
The (absolute value of the) charge of each particle is 
x/qfes, the square root of the Fermi edge singular- 
ity exponent of the model. This quantity is defined 
through the behavior of the zero-temperature Green 
function of d^ip(0) for 7T W = 0, which, for long time 
t decays as T ~ aFES . We will give expressions for this 
quantity below. For a wire of a finite length L at zero 
temperature, the Coulomb gas interaction is given by 
In | sinh [J!ify{_^..)/x,] | • At finite temperatures it takes 
the form of a logarithm of an elliptic function [1]. The 
other part of the Coulomb gas Hamiltonian is simply 
equal to the total length of the imaginary time inter- 
vals in which the level is occupied, multiplied by eo, 
the energy of an occupied level. It is analogous to an 
electric field applied on the classical system of charges. 

Up to now, the values of «fes and r have not been 
specified. When there are no intra-wire interactions 
(i.e., g — 1), we have the usual interacting resonant 
level model, in which case qfes = (l — §5) , and 
To = it |7i w | pocos(5), where 8 = tan _1 (7rpoWi w /2) is 
the phase shift experienced by electrons at the Fermi 
surface in the wire due to the level-wire interaction, 
and po is the local density of states in the wire edge at 
the Fermi energy [17,18]. For an interacting wire (g ^ 
1) the situation is somewhat more complicated. Since 
there is no backscattering in this problem, calculations 
using bosonization [1] yield Qfes = (1 — gU\ w /irv a ) 2 / g, 
and To = 7r|7i w | 2 po- Comparing these expressions 
with the previous ones in the limit of g — 1 (noninter- 
acting wire), we see that the results of bosonization 
replace the phase shift 5 by its first Born approxima- 
tion value. Thus, the values of «fes and To in any par- 



ticular systems are renormalized by irrelevant terms 
not appearing in the Tomonaga-Luttinger Hamilto- 
nian (1). It is natural to expect that, taking these 
effects into account, for a general model we would 
get qfes = j (l " <^eft) , for some effective phase 
shift 5 e s G [— it/2, it/2], so that F will be given by 
7T |7i w | 2 vo cos(5 e ff)- This can be shown to be valid for 
models in which ofes can be exactly evaluated [13]. 
Finally, it should be noted that in general eo is also 
modified by a term representing the difference in the 
total energy of the wire caused by the potential applied 
on it by an empty or a filled level. This correction can 
also be related to the phase shifts. Its value, however, 
vanishes at half filling due to particle-hole symmetry, 
which is the case in our numerical calculations. 

Let us now turn to the LDoS, which we shall denote 
by pievei(w). This quantity is equal (up to a factor of 
— 1/ir) to the imaginary part of the level retarded Green 
function, which in turn can be obtained by analytic 
continuation from the corresponding Matsubara Green 
function in the upper complex frequency plane [19]. In 
the imaginary time domain the Matsubara Green func- 
tion is defined by Gi(r) = -Tr{f T e- H/T a(r)a t (0)}/Z, 
where T T is the imaginary time ordering operator. The 
numerator of this expression can be given a Coulomb- 
gas representation, which has the same form as Eqs. (3- 
4) with two additional charges of sizes ±(^ofes — 
1/y/g), inserted at r and at the origin, respectively. 
These charges correspond to the level creation and an- 
nihilation operators appearing in the definition of the 
Green function. 

From the above results we can immediately see an 
interesting distinction between the Coulomb-gas ex- 
pression for the partition function and that for the 
LDoS: The former contains only three parameters: To, 
eo, and Qfes, while the latter explicitly depends on g 
too. As both the interactions in the wire and the level- 
wire interactions affect the partition function mainly 
though a single combination — the Fermi edge singu- 
larity exponent oifes, thermodynamics cannot be used 
to distinguish between the different interaction types. 
In other words, one can construct very different mod- 
els, whose interactions differ in strength and even in 
sign, which will have the same thermodynamic prop- 
erties, provided To, eo, and qfes are indeed the same. 
On the other hand the LDoS, which depends explicitly 
on g, can be used to disentangle the effects of intra- 
wire and level-wire interaction, and will thus behave 
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Non-interacting lead, U |w ^0 
Interacting lead, U |w =0 
Interacting lead, U. ^0 




Fig. 1. (Color online) DMRG results for the level popula- 
tion as a function of its energy, for three different models 
denoted by the three different symbol types. The curves on 
which the symbols reside (which serve as a guide to the eye) 
correspond to the various (Jfes values (the larger ofes the 
narrower the curve and vice versa) . See the text for further 
details. 

differently for these different systems. In fact, it can 
be expected that at low energies the LDoS should be- 
have like the tunneling density of states near the end 
ofaLL, i.e., it will vary as [max(|w|, T, v/L)] 1/s ' L [1]. 
Thus, it exhibits a power law, which depends only on 
the LL parameter g and not on qfes, i.e., on the inter- 
actions in the wire but not on the level-wire coupling. 
These considerations thus give rise to a quite surprising- 
possibility: In spite of the fact that the level popula- 
tion is the integral of the LDoS times the Fermi-Dirac 
distribution function, a LL system (i.e., a system with 
g / 1) may have the same dependence of the level oc- 
cupancy on £o as a Fermi liquid system (for which g = 
1) with a properly chosen level-wire coupling, although 
the LDoS is obviously different in the two cases. In the 
next section we will demonstrate this to hold quanti- 
tatively, using numerical calculations. 



3. Numerical results 

In this section we show numerical data confirming 
the results of the previous section, i.e., that the level 
population is universal (equal for different models), al- 
though the LDoS is not. 

To find the level occupation we have performed den- 
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Fig. 2. (Color online) Monte Carlo results for the LDoS as a 
function of frequency, measured from the Fermi energy for 
three different models, all with ctpES = 2/3, eo/Fo = 0.7, 
and T/To = 0.04. The three curves correspond to three dif- 
ferent values of g, as indicated in the legend. Inset: temper- 
ature dependence of the LDoS at the Fermi energy (sym- 
bols), together with the expected power-low behavior (line). 
See the text for further details. 



sity matrix renormalization group (DMRG) [20] calcu- 
lations on a particular realization of a LL wire repre- 
sented by a half-filled iV-site tight binding chain with 
nearest neighbor interactions. In this case, the wire and 
level-wire Hamiltonians are given by: 



Hw = (— * c I c i+i +H-c.) + 



U ( c]a - t 



d J d-- 



,t 

-i+i 



Ci+i - g 



#iw = - (ti v ,c[d + H.c.J + 
Ui 



t 1 

4d - 2 



(5) 



(6) 



where c\ (a) is a creation (annihilation) operator for an 
electron at the wire's i'th site, t and U are the nearest- 
neighbor hopping and interaction strengths along the 
chain, and ti w , Ui w denote the corresponding quantities 
in the level-wire coupling term. The latter are related 
to the parameters of the continuum level-wire coupling 
Hamiltonian (2) by 7T W = t\ w y/a, and ti\ w = U\ w a, a 
being the lattice spacing. Using boundary conformal 
field theory arguments and the Bethe ansatz, it can be 
shown that for this model [13] : 
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Qfes = - 
9 



i 2 5 i 

1 - tan , 

7T \ y/(2t) 2 - U 2 



V 



(7) 



where g = tt/[2 cos _1 (-I//2t)]. 

In Fig. 1 we show the level population as a function 
of its energy. The different curves correspond to differ- 
ent «fes values. On each such curve there are symbols 
of three types, denoting DMRG data in three different 
systems: (i) U = but f7i w / 0; (ii) U ± but (7i w = 0; 
(iii) both U ^ and {7i w / 0. The interactions in the 
three three models were chosen so as to give equal ofes 
values, as denoted in the figure's legend [For model (c) 
we used U = ±0.5t, with sign opposite to that of model 
(b)]. In all cases we have kept To = 10~ 4 i (by choosing 
ti w appropriately), and N = WOv/t. A block size of 
256 was used. From the results one can clearly see that 
the population is indeed universal, and is equal for dif- 
ferent models provided oifes is the same, although the 
interaction strengths are different in magnitude and 
sign for the three cases. 

From the data of Fig. 1 we observe that the pop- 
ulation vs. level energy curves become wider as «fes 
becomes smaller and vice-versa. Our previous results 
show that smaller ofes implies either g > I (i.e., at- 
tractive interactions in the wire) or J7i w > 0. This could 
be understood since both options should enhance the 
effective level-wire hopping, and thus lead to a wider 
population curve. Indeed, for g > 1 it is well known 
that the local density of states at the edge of a LL 
(or at the middle of a chiral LL) is enhanced at the 
Fermi energy [1], so level- wire hopping becomes effec- 
tively stronger. Similarly, f7i w > also facilitates larger 
tunneling by the Mahan exciton effect [17]: due to the 
level-wire repulsion, when the level is empty the site 
at the wire's end tends to be occupied and vice versa; 
hence, it helps the electrons overcome the limitations 
on hopping induced by the Pauli principle. In the same 
way we can understand why larger ofes will cause nar- 
rower population curves. 

To find the LDoS we used classical Monte-Carlo 
calculations to compute the imaginary-time Green 
function from its Coulomb gas representation. The 
results were then Fourier transformed into Matsub- 
ara frequency domain, and analytically continued to 
obtain the retarded Green function using the Pade 
approximant technique [21]. The results are shown in 
Fig. 2. The three curves have parameters which are 
approximately equal to those of the three models rep- 



resented by the widest curve of Fig. 1. We immediately 
see that, although the populations are, to a very high 
accuracy, equal in the three models, since ofes = 2/3 
for all of them (actually, the Coulomb gas representa- 
tion would predict exactly identical occupations), the 
LDoS are markedly different, having a maximum, a 
minimum, or no special feature near the Fermi energy 
for g > 1, g < 1, and g = 1, respectively. In the inset 
we demonstrate that the LDoS at the Fermi energy 
has a power-law dependence on temperature for all 
three cases, with the expected power of 1/g — 1. 



4. Conclusions 

To conclude, we have shown, both analytically and 
numerically, that the thermodynamics of a level cou- 
pled to the edge of a LL is universal, depending on only 
a few parameters. Thus, thermodynamic quantities (in 
particular — the level occupation) of a level coupled 
to a LL can be equal to the corresponding quantities 
in other systems which are Fermi-liquids. This occurs 
although thermodynamic quantities are determined by 
dynamic properties which are not universal. In this 
work we have concentrated on the level occupation, 
which is determined by the LDoS. The latter exhibits 
a power law behavior at low frequency with a power 
directly related to LL parameters, which cannot be re- 
produced by a Fermi- liquid. Yet, this LL-specific power 
law does not reflect itself in a LL-specific ni cvo i(eo) de- 
pendence. Thus, thermodynamics of quantum impuri- 
ties may not help to expose LL physics. Nevertheless, 
our results imply that interesting phenomena in such 
systems may be studied on equivalent models with non- 
interacting wires, which are much easier to approach, 
both analytically and numerically (using, e.g., Wilson's 
numerical renormalization group [22]). 
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